Естественные науки 


УДК 537.874.4 


ИССЛЕДОВАНИЕ ЭЛЕКТРОМАГНИТНОГО РАССЕЯНИЯ СТРУКТУРАМИ, 
СОСТАВЛЕННЫМИ ИЗ НЕСКОЛЬКИХ НЕПЕРЕСЕКАЮЩИХСЯ ПРОВОДНИКОВ 

Ю.А. Келлер 

Томский государственный университет 
Е-плаіІ: киа1102@1ртк.1$и.ги 

На основе метода вспомогательных источников построен численный алгоритм решения задач электромагнитного рассеяния на 
структурах, составленных из конечного числа непересекающихся тонких проводников. Построенный алгоритм реализован в ви- 
де компьютерной программы для расчета характеристик рассеяния ряда структур, отличающихся взаимным расположением 
проводников. Исследовано влияние взаимного расположения проводников на бистатические сечения рассеяния рассмотренных 
структур, а также на распределения тока вдоль проводников. 


Теория возбуждения и рассеяния электромаг- 
нитных волн тонким проводником родилась давно 
[1]. На начальном этапе эта теория была порождена 
запросами антенной техники. Однако, начиная с се- 
редины 60-х гг. прошлого века, исследователи начи- 
нают уделять всё большее внимание анализу рассе- 
ивающих свойств тонких проводников. Интерес к 
подобным исследованиям был обусловлен потреб- 
ностями создания объектов с заданными рассе- 
ивающими свойствами. Наиболее распространен- 
ным методом решения как задач возбуждения тон- 
кого проводника, так и задач рассеяния на нём, яв- 
ляется метод интегральных уравнений (см., напри- 
мер, [2-4]). Если рассеивающая структура состоит 
из нескольких проводников, то решение рассматри- 
ваемых задач сводится к решению системы ней- 
тральных уравнений, порядок которой равен числу 
проводников в структуре. В данной работе для ре- 
шения задач рассеяния на структурах, содержащих 


тонкие проводники, использован вариант метода 
вспомогательных источников, который позволяет 
исключить этап построения системы интегральных 
уравнений, тем самым упрощая процесс решения 
задачи и снижая затраты в использовании компью- 
терных ресурсов. Даны математическая формули- 
ровка варианта и краткое описание возможностей 
реализованной на его основе программы для расчё- 
та распределений тока вдоль проводников и харак- 
теристик рассеянного поля различных структур. 
Приведены результаты численных результатов, ха- 
рактеризующих взаимное влияние проводников на 
распределение тока вдоль них и бистатические сече- 
ния рассеяния (БСР). 

1. Формулировка задачи и метод её решения 

Геометрия задачи показана на рис. 1. Будем рас- 
сматривать стационарную (зависимость от време- 



Рис. 1 . Г еометрия задачи 
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ни выбрана в виде ехр (—Ш\ задачу дифракции 
электромагнитного поля { Д, Щ на структуре, со- 
стоящей из II непересекающихся проводников, 
ограниченных поверхностями Д (и= 1 ,2,..., II) и рас- 
положенных произвольным образом по отноше- 
нию друг к другу. Под тонким проводником будем 
понимать идеальный проводник круглого сечения, 
диаметр которого мал по сравнению с длиной вол- 
ны и длиной проводника. Эта структура размещена 
в однородной безграничной среде Д. с диэлектри- 
ческой и магнитной проницаемостями в е и /и е соот- 
ветственно в декартовой системе координат с нача- 
лом отсчета, лежащем внутри проводника с поряд- 
ковым_номером и= 1. Требуется найти рассеянное 
поле {Е е ,Щ в области Д. 

Математическая постановка задачи имеет сле- 
дующий вид: 


[V, Д ] = ісоц е Н е , [Ѵ,Н е ]=ісое е Е е 

(1) 

в области Д, 


[Д Д] = -[»„, Д] 

(2) 

на поверхностях Д, где и=1,2,..., Д 


[^Ё'Л/К] +4Кн е =0{к-\ 


[^ е Н е Я/К]-^Ё е =0{К- 1 ) 

(3) 


при Л— > со. Здесь л„ - единичные вект оры нор малей 
к поверхностям Д проводников, Л=л/хЧуЧД. 

Решение сформулированной выше задачи по- 
лучим следующим образом. Разместим (рис. 1) 
внутри каждого из тонких проводников на его оси 
непрерывно распределенный вспомогательный ток 
7 и . Представим неизвестное рассеянное поле {Е„Щ 
в области Д в виде суммы полей введенных вспо- 
могательных токов: 

Дд^) = й[ѵ,[ѵ,п„]], 

к е и = 1 

н е (М) = — Х[Ѵ,П„], 

Ие М =1 

П а =^ е (М,М Іл йЛ (4) 

І и 

Здесь к=со7ё]7 г , Ѵ е (ММіи)=(4яК т )-'ехр(ікАщ)', 
К М Мім — расстояние от точек М І и на оси проводников 
до точки наблюдения М в области Д; 7 и - неизвест- 
ные осевые вспомогательные токи, и=1,2,..., Д инте- 
грирование проводится вдоль осевых линий провод- 
ников; /„ - осевая линия проводника с номером и. 

Поле (4) удовлетворяет уравнениям Максвелла 
( 1) и условиям излучения (3) в области Д. Для того, 
чтобы удовлетворить граничному условию (2), 
необходимо соответствующим образом выбрать 
осевые вспомогательные токи 7 и (и= 1,2,..., Ц). 

Введем кусочно-постоянную аппроксимацию 
осевых токов. Разобьем линию /„ каждого тока 7„ на 
/Ѵ„ малых участков, в пределах каждого из которых 
ток можно считать постоянным. Тогда выражение 
для П„ в (4) приближенно можно записать в виде 


П„ = ІЛЛ, }Ч(Л/,ЛД„)Д (5) 

М І.М 

где Д — ток на /-ом участке проводника с номером 
и, е иі - единичный вектор, направление которого 
совпадает с направлением касательной в средней 
точке рассматриваемого участка. При таком подходе 
нахождение неизвестных распределений осевых то- 

V 

ков сводится к нахождению элементов тока. 

и = 1 

Для определений значений элементов тока ис- 
пользуем граничные условия (2), удовлетворяя им в 
соответствии с методом коллокаций. Пусть М і 
(/=1,2,...,Д) - точки коллокации на поверхности 
Д; Д - число точек коллокации на Д. В силу пред- 
положения о малости диаметра проводника по 
сравнению с длиной проводника и длиной волны 
будем считать, что вкладом в рассеянное поле ази- 
мутальных составляющих токов на поверхностях 
тонких проводников можно пренебречь. Тогда для 
нахождения неизвестных элементов токов Д 
(и=1,2,...,Д /=1,2,...Д,) получим следующую си- 
стему линейных алгебраических уравнений с ком- 
плексной матрицей размерности Д ) х /Ѵ ц ) : 

и= 1 и= 1 

Ки>=-Ки» «=и-..ДЛ У = 1,2,..., Д, (6) 

где Д„, и - значения электрических компо- 
нент рассеянного (4) и возбуждающего полей вдоль 
оси проводника с номером и в точках коллокации 
на его поверхности. 

Решение системы (6) определяется путём мини- 
мизации функционала 

ф =ііѢ Е и 1 +к, І і 2 - < 7 > 

М-1 7=1 

После решения задачи минимизации (опреде- 
ления неизвестных элементов тока Д , (и=1,2,..., II; 
і=1,2,...,Щ необходимые характеристики рассеян- 
ного поля определяем из (4). В частности, для ком- 
понент рассеянного поля в дальней зоне имеем 

Е еіѲ (М) = ( Ие /8 е У /2 Н е№ (М) = 

= (ехр(ДД)/ ДК)Д(0,<р) + 0(К Д, 

Е С9 (М) = -(^/е с )' /г Н ер (М) = 

= (ехр (ік'К)/к'№ 9 (Ѳ,<р) + 0(К 2 ), (8) 

где К - расстояние от начала системы координат до 
точки наблюдения М, а Д(0,(р), Д(0,<р) - компоненты 
диаграммы рассеяния, определяемые выражениями 

и N и 

щѳ,Ф)= 2ХЛ ДС08ѲС08(рС08а м/ . + 

и= 1 /=1 

+ С08 Ѳ 8ІП ср С08 р иі - 8ІП Ѳ С08 у и . ) І ц і , 

V Ы и 

°УѲ,(р) = ^^7 иі (-8ІпѲсо8а и . + 

и= 1 /=1 

+со8<рсо$/Д,.)Д,, (9) 
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где 7„,=і<у/„, - решение си- 

стемы (6), со8а иі , со8/3„„ со$у„, - направляющие 
косинусы единичного вектора е иі , Ѳ и (р - общепри- 
нятые угловые сферические координаты точки на- 
блюдения М, а интеграл І иі имеет следующий вид: 

*А,і 

I иі = [ ехр{-/( ып 0 сов <рк е х, и + 

+ &тѲ$іп(рк е у Іи + со &Ѳк е і Іи )}сІ(к е І), 

где крСі м , к^ іи , к е ъ і<и — координаты на оси проводни- 
ка, по которым проводится интегрирование. 

Данный метод, как и другие варианты метода 
вспомогательных источников, позволяет осущест- 
влять апостериорную оценку точности полученно- 
го решения. В качестве величины, характеризую- 
щей точность, выберем значение относительной 
нормы невязки граничных условий на поверхно- 
стях всех проводников 5 и в точках, промежуточных 
по отношению к точкам коллокации 

л = (ф'/ф 0 ) 1/2 > ф„ = ІІі</і 2 > О») 

и = 1 ] =1 

где Ф ’ - значение функционала (7) на указанной 
выше совокупности точек; Ф 0 - значение нормы 
падающего поля на этой же совокупности точек; 
Ь' и - число промежуточных точек на поверхности 
проводника с номером и. 

2. Численные результаты 

На основании изложенного выше метода созда- 
на компьютерная программа для расчета компо- 
нент рассеянного поля и контроля точности полу- 
ченного решения для структур, конфигурация ко- 
торых представлена на рис. 2-4. 

Первая структура (рис. 2) представляет собой 
одиночный прямолинейный проводник длиной Ь ь 
расположенный в декартовой системе координат 
ОХУ 2 таким образом, что его осевая линия напра- 
влена вдоль оси 02, а середина осевой линии сов- 
падает с началом системы координат. Вторая 
структура (рис. 3) состоит из двух параллельных 
проводников длиной Ц и Ь 2 с осевыми линиями, 
направленными вдоль оси 02, расположенных на 
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расстоянии 5 друг от друга; середина осевой линии 
второго проводника расположена на оси ОХ. 
Третья структура (рис. 4) состоит из взаимно орто- 
гональных проводников. Осевая линия первого 
проводника этой структуры длиной X, ориентиро- 
вана вдоль оси 02, её середина совпадает с нача- 
лом системы координат. Два других проводника 
длиной 2, и 2 3 расположены на одинаковых рас- 
стояниях <5 от первого проводника таким образом, 
что их осевые линии направлены вдоль оси ОУ. 

Входными величинами программы являются 
конфигурация структуры, поле {Д,Д}, длины и ра- 
диусы проводников, расстояние между проводни- 
ками, а также число элементов разбиения осевого 
тока /Ѵ„ для каждого из проводников. 

Минимизация функционала (7) осуществляется 
методом сопряженных градиентов; итерационный 
процесс останавливается при условии, если изме- 
нение функционала после выполнения очередной 
итерации не превышает КУ 5 . 

При помощи данной программы выполнена се- 
рия вычислений, направленных на выяснение 
влияния числа элементов разбиения осевого тока 
на величину невязки граничных условий, на срав- 
нение получаемых результатов с результатами дру- 
гих авторов, а также на оценку взаимного влияния 
проводников на распределения тока вдоль них и 
БСР составленных из них структур. Некоторые ре- 
зультаты представлены ниже. 

Предполагается, что структуры возбуждаются 
таким образом, что векторы Е п и к е лежат в плоско- 
сти 20У и вектор к е образует с осью ОУ угол у/ (как 
показано на рис. 2-4). Предполагается также, что 
радиус проводников во всех случаях равен 0,02А, 
где А - длина возбуждающей волны. 

На рис. 5 представлены результаты сравнения 
распределения осевого тока вдоль одиночного про- 
водника длиной 2= А при угле падения плоской 
волны у/=30°. 


И, уел. ед. 



Рис. 5 . Распределение осевого тока вдоль проводника дли- 
ной 2= А 


По оси абсцисс отложена координата точки на 
оси проводника, а по оси ординат - значение тока в 
этой точке в условных единицах. Кривая 1 - резуль- 
таты, полученные данным методом, кривая 2 - ре- 


зультаты из работы [5] (получены методом интеграль- 
ных уравнений). При получении кривой 1 линия то- 
ка разбита на 40 участков. При данном числе разбие- 
ний значение невязки равно 0,284. Как показывают 
эти кривые, имеет место хорошее совпадение срав- 
ниваемых результатов. Незначительные наблюдае- 
мые отличия можно объяснить ошибками при гра- 
фическом съёме информации с рисунка работы [5], а 
также погрешностями вычислений как данным мето- 
дом, так и методом интегральных уравнений. 

Результаты, представленные на рис. 6-9, позво- 
ляют оценить взаимное влияние параллельных 
проводников на токовые распределения и БСР. Ис- 
следования такого рода представляют интерес для 
теории и техники проволочных (вибраторных) ан- 
тенн, а также для оценки радиолокационной за- 
метности объектов, содержащих такие провода. 
Оценка взаимного влияния проводников осущест- 
влялась путём сравнения токов распределений 
вдоль проводников и БСР структур, состоящих из 
двух (рис. 3) параллельных проводников, с соответ- 
ствующими характеристиками для одиночного 
проводника (рис. 2). 

Рис. 6 и 7 характеризуют распределение токов 
вдоль оси проводников длиной 2=А, расположен- 
ных очень близко (<5=0,0016А) и достаточно далеко 
(<5=ЗА) друг от друга. 


|7|,усл.ед. 



Рис. 6. Распределение токов вдоль оси проводников длиной 
2= А, расположенных на расстоянии <5=0,00 16А 


ІДусл.ед. 



Рис. 7 . Распределение токов вдоль оси проводников длиной 
2= А, расположенных на расстоянии 8= ЗА 
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Проводники возбуждаются плоской волной, 
падающей под углом і//=30°. По оси абсцисс отло- 
жена координата точки на оси проводника, по оси 
ординат - значение тока в этой точке в условных 
единицах. Значение к е 1= 0 соответствует середине 
проводников, значения к е [> 0 относятся к верхней 
части проводников, значение к е І<0 - к нижней ча- 
сти проводников. Кривые 1 на рис. 6, 7 - распреде- 
ление тока вдоль одиночного проводника; кри- 
вые 2 - распределение тока вдоль каждого из про- 
водников системы из двух проводников (распреде- 
ления одинаковы). Во всех случаях число элемен- 
тов разбиения осевого тока выбиралось равным 40. 

На рис. 8, 9 представлены БСР 
а(Ѳ,<р) = 1іш{| Е ев {Ѳ,ср) | 2 + 1 Е„(Ѳ іЧ >) | 2 }/ 1 Ё 0 \ 2 } (11) 

К-*со ’ 


тех же структур в полуплоскости ср=0° для взаимно- 
го расположения проводников: <5=0,0016А (рис. 8) 
и <5=ЗА (рис. 9). 



Рис. 8. БСР структур, состоящих из проводников длиной 
Ь= А, расположенных на расстоянии <5=0,0016А 


7ж дБ - 



Рис. 9. БСР структур, состоящих из проводников длиной 
І=А, расположенных на расстоянии <5=ЗА 


Структуры возбуждаются плоской волной, па- 
дающей под углом і//=30°. На этих рисунках кривые 
1 - БСР одиночного проводника (рис. 2), кривые 2 
- БСР структуры из двух параллельных проводни- 
ков (рис. 3). 


Представленные на рис. 6-9 результаты расчё- 
тов позволяют сделать следующие выводы. 

Если вблизи прямолинейного проводника на- 
ходятся другие проводники, то распределение то- 
ков вдоль проводников существенно отличаются от 
распределения тока вдоль одиночного проводника 
(рис. 6). Если структура состоит из двух проводни- 
ков, то распределение тока на обоих проводниках 
одинаковы. Отличие токовых распределений для 
структур проводников при расположении их близ- 
ко друг к другу от распределения тока вдоль оди- 
ночного проводника объясняется сильным взаимо- 
действием этих проводников. Однако, несмотря на 
существенные различия токовых распределений 
структур близко расположенных проводников от 
распределения тока вдоль одиночного проводника, 
БСР одиночного проводника и структур из двух 
проводников отличаются мало (рис. 8). Это гово- 
рит о том, что по отношению к рассеянному полю 
структура из близко расположенных параллельных 
проводников эквивалентна одному проводнику. 
Если проводники рассматриваемых структур нахо- 
дятся достаточно далеко друг от друга (в данном 
случае на расстоянии <5=ЗА), то распределения то- 
ков вдоль проводников близки к распределению 
тока вдоль одиночного проводника (рис. 7). Это 
объясняется значительным уменьшением взаимо- 
действия проводников структуры. В этом случае 
можно считать каждый проводник независимым 
рассеивателем, и рассеянное поле, следовательно, 
является суперпозицией полей, рассеянных каж- 
дым проводником, т.е. имеет интерференционную 
структуру, что и отражают результаты, рис. 9. 

Результаты, представленные на рис. 10, 1 1, позво- 
ляют оценить взаимное влияние ортогональных про- 
водников на токовые распределения и БСР. Как и в 
случае параллельных проводников, оценка взаимно- 
го влияния осуществлялась путём сравнения токо- 
вых распределений и БСР структуры, состоящей из 
перпендикулярных проводников (рис. 4) с соответ- 
ствующими характеристиками для одиночного про- 
водника (рис. 2). При проведении численных расче- 
тов во всех случаях число элементов разбиения осе- 
вого тока выбиралось равным 40 на длину волны. 



Рис. 10. Распределение тока вдоль первого проводника струк- 
туры при падении плоской волны под углом і//=30° 
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Рис. 11. БСР структуры из взаимно перпендикулярных провод- 
ников при падении плоской волны под углом і//=30° 

На рис. 10 показаны распределения токов вдоль 
первого проводника структуры, представленной на 
рис. 4 при падении плоской волны под і//=30° и раз- 
личных расстояниях между проводниками. Длина 
всех проводников структуры равна Я/2. Значение 
к е 1= 0 соответствует середине проводников; значе- 
ния к е [> 0 относятся к верхней части вертикального 
(первого) проводника структуры; значения к е 1<0 - к 
нижней части вертикального проводника. 

Кривая 1 характеризует распределение тока 
вдоль вертикального (первого) проводника струк- 
туры при (5=0, 05Я, кривая 2 - при <5= 0,2Я, кривая 3 
- при <5=0, 8Я и кривая 4 - при <5=1, 5Я, кривая 5 ха- 
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рактеризует распределение тока вдоль одиночного 
вертикального проводника. 

На рис. 11 изображены БСР в полуплоскости 
<р=0° рассматриваемой структуры для того же угла 
падения плоской волны и расстояний между провод- 
никами. Обозначения кривых аналогичны рис. 10. 

Выводы 

На основе метода вспомогательных источников 
построен численный алгоритм и реализована ком- 
пьютерная программа решения задач электромаг- 
нитного рассеяния на структурах, составленных из 
конечного числа непересекающихся тонких про- 
водников. Исследовано влияние взаимного распо- 
ложения проводников на бистатические сечения 
рассеяния рассмотренных структур, а также на ра- 
спределения токов вдоль проводников. 

При наклонном падении волны к оси централь- 
ного проводника структуры падающая волна возбуж- 
дает как центральный проводник, так и расположен- 
ные рядом перпендикулярно к нему боковые провод- 
ники. В этом случае при малых расстояниях между 
проводниками (<5<0,2А) токовые распределения на 
центральном проводнике зависят от расстояния; од- 
нако при расстояниях <5>0,8А токовые распределения 
на центральном проводнике структуры мало отлича- 
ются от распределения тока вдоль такого же одиноч- 
ного проводника. Показано, что бистатические сече- 
ния рассеяния рассматриваемой структуры отлича- 
ются от таковых для одиночного проводника. 
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